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Abstract 

It is shown that the ratio of strange to non-strange particle production in rela- 
tivistic heavy ion colhsions is expected to reach a maximum at beam energies around 
30 A-GeV in the lab frame. This maximum is unique to heavy ion collisions, and 
has no equivalent in elementary particle collisions. The appearance of the maxi- 
mum is due to the energy dependence of the chemical freeze-out parameters and 
is clearly seen as a pronounced peak in the Wroblewski factor as a function of the 
incident energy as well as in the behavior of the /tx^ ratio. Below 30 A-GeV 
strange baryons contribute strongly because of the very large value of the chemical 
potential ^b- As the energy increases, the decrease of the baryon chemical potential 
coupled with only moderate increases in the associated temperature causes a decline 
in the relative number of strange baryons above energies of about 30 A-GeV leading 
to very pronounced maxima in the A/tt^ and S"/??"*" ratios. 



1 Introduction 



The experimental data from heavy ion colhsions show that the j-n^ ratio is larger at BNL- 
AGS energies than at the highest CERN-SPS energies [1-6] and even at RHIC [0. This behavior 
is of particular interest as it could signal the appearance of new dynamics for strangeness 
production in high energy collisions. It was even conjectured that this property could indicate 
an energy threshold for quark-gluon plasma formation in relativistic heavy ion collisions. In 
this paper we analyze the energy dependence of strange to non-strange particle ratios in the 
framework of a hadronic statistical model. The statistical approach has been very successful in 
describing particle yields from low energies starting with SIS and AGS all the way up to SPS 
[9-16] and RHIC [^ |T^, |T^ energies. In particular it was found that, in the whole energy range, 
the hadronic yields observed in heavy ion collisions resemble those of a population in chemical 
equilibrium along a unified freeze-out curve determined by the condition of fixed energy/particle 
~ 1 GeV [|^. The chemical freeze-out curve provides a relation between the temperature T 
and the baryon chemical potential /x^ when the particle composition of the system is frozen in. 

In this paper we use the energy dependence of the thermal parameters at chemical freeze- 
out as determined from the analysis of [10-16] to show that the ratio of strange to non-strange 
particle multiplicities reaches a maximum at approximately 30 A- GeV lab energy. Compared to 
other qq pairs, the relative strangeness content actually declines by almost a third between 30 
A- GeV and RHIC energies. We emphasize that this maximum is unique to heavy ion collisions 
and does not occur in the collisions of elementary particles like p-p, p-p or e"'"-e~. 

The values of the baryon chemical potentials and temperatures at different collision energies 
are shown in Fig. 1. The values of /i^ at SIS, AGS and SPS energies were taken from a detailed 
analysis [|I^, 13] of the data on particle multiplicities extrapolated to full phase space while the 



value of /is at RHIC was obtained from data taken at mid-rapidity |jT8[. The energy dependence 
can be parametrized phenomenologically as 

where a ~ 1.27 GeV and h ~ 4.3 GeV. The results of this parametrization are shown by the 
full line in the upper part of Fig. 1. At freeze-out the chemical potential is related to the 
temperature via the phenomenological condition of fixed energy per hadron [1T3, namely 



(E) / (iV) ^ 1 GeV. (2) 

Applying Eq. (1) together with Eq. (2) leads to the energy dependence of the temperature 
shown by the full line in lower part of Fig. 1 together with the values obtained in [|T^, ^ . 

Both the temperature and the chemical potential exhibit a strong variation with energy. 
The GSI/SIS results have the lowest freeze-out temperatures and the highest baryon chemical 
potentials. As the beam energy is increased a clear shift towards higher T and lower /i^ occurs. 
Above AGS energies, the temperature exhibits only a moderate change and converges to its 
maximal value in the range of 160 to 180 MeV. This value is remarkably close to the critical 



temperature Tc ~ 170 ± 8 MeV extracted from the lattice calculations of QCD ||19 . 

In the statistical approach, the energy dependence of the basic thermal parameters deter- 
mines the energy dependence of many relevant observables, for instance particle yields. For 
strange particle production this determines the ratios of strange to non-strange particle multi- 



plicities as well as the Wroblewski factor [^ defined as 



' ~ {uu) + {dd) ^ ^ 
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where the quantities in angular brackets refer to the number of newly formed quark-antiquark 
pairs, i.e., it excludes all quarks that were present in the target and projectile. 

Applying the statistical model to particle production in heavy ion collisions calls for the use 
of the canonical ensemble to treat the number of strange particles particularly for data in the 
energy range from SIS up to AGS |T^, ^ . For these energies, the number of strange particles per 
event is so small that the description of strangeness conservation on the average, as contained 
in the grand canonical approach is not adequate. The exact conservation of quantum numbers 
in relativistic statistical mechanics has been well established for some time now ||22|, E3| . In the 



following we present a new analytical expression of the canonical partition function and the 
corresponding multi-strange particle multiplicities. This is quantitatively equivalent to previous 
results |2^, however, the result can now be presented in a compact format well suited for 
numerical evaluation. 



2 Canonical strange particle multiplicities 

The canonical partition function of a hadronic resonance gas constrained by the strangeness 



neutrality condition reads P3 
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exp Sne''^'^ ' 



(4) 



where Sn = V X]fc Z], with V being the volume and the sum is over all particles and resonances 
carrying strangeness n. For a particle of mass m^, with spin-isospin degeneracy factor gk, 
carrying baryon number Bk and charge Qk-, with the corresponding chemical potentials and 
/xq, the one-particle partition function is given, in Boltzmann approximation, by 
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The integral representation of the partition function in Eq. (4) is not convenient for a 
numerical analysis as the integrand is a strongly oscillating function. However, the (p integration 
in Eq. (4) can be done exactly. Indeed, rewriting Eq. (4) as 
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and using the relation 6 2^*+*) = Y.-'^t'^Imip) one obtains, after integrating, the partition 
function in a form free of oscillating terms 



where 
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(7) 



(8) 
(9) 



and li are modified Bessel functions. The contribution of thermal pions to 5*0 is calculated 
using Bose-Einstein statistics. 
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The expression for the particle density, nj, can be obtained from the partition function 
Eq. (3) by following the standard method pSf. For a particle i having strangeness s the result 
is 

^1 oo oo 

Hi = ala2a'^'^''~^^~'In{x2)Ip{x^)I-2n-zp-s{xi). (10) 

5=0 n=— oop=— oo 



The leading term in X3 corresponds to the approximate form used in ||2^ to study the centrality 
dependence of strangeness production in Pb-Pb collisions at SPS energies. 

It can be verified that, in the limit of large Xj, the above formulae coincide with the grand 
canonical result. In the opposite limit, that of small Xj, the equilibrium density of (multi) strange 
particles is strongly suppressed relative to its grand canonical value. The above equation can 
therefore be applied both at SIS energies where the number of strange particles produced per 
event is as small as 10~^, and at SPS where over a hundred strange particles are produced on 
average in each event. 



3 Results 

The theoretical results of the previous section combined with the energy dependence of the 
thermal parameters shown in Fig. 1 provide a basis for the study of the energy dependence of 
strangeness production in heavy ion collisions. Of particular interest are the ratios of strange 
to non-strange particle multiplicities as well as the relative strangeness content of the system 



as expressed by the Wroblewski factor ||20|. The thermal model calculations for Au-Au and 
Pb-Pb collisions are performed using a canonical correlation volume given by a radius of ~ 7 
fm. This radius could vary with energy but for simplicity this was not taken into account 
in our calculations. At a later stage, when more data become available, this approximation 
will of course have to be reconsidered. Furthermore, we assume strangeness to be in complete 
equilibrium, that is, the strangeness saturation factor is taken as 7s = 1. When more data will 
become available, this point might be refined. 

We turn our attention first to the energy dependence of the Wroblewski ratio The quark 
content used in this ratio is determined at the moment of chemical freeze-out, i.e. from the 
hadrons and especially, hadronic resonances, before they decay. This ratio is thus not an easily 
measurable observable unless one can reconstruct all resonances from the final-state particles. 
The results are shown in Fig. 2 as a function of invariant energy ^/s. The values calculated 
from the experimental data at chemical freeze-out in central A-A collisions have been taken 



from reference [|T5[[]. The values of were extracted from fully integrated 47r data with the 



exception of the result form RHIC where particle ratios measured at mid-rapidity [Q were used. 
The solid line in Fig. 2 describes the statistical model calculations in complete equilibrium along 



the unified freeze-out curve |jT^ and with the energy dependent thermal parameters presented 
here. From Fig. 2 we conclude that around 30 A-GeV lab energy the relative strangeness 
content in heavy ion collisions reaches a clear and well pronounced maximum. The Wroblewski 
factor decreases towards higher incident energies and reaches a limiting value of about 0.43. 

The appearance of the maximum can be traced to the specific dependence of on the 
beam energy. In Fig. 2 we also show the results for calculated under the assumption that 
only the temperature varies with collision energy but the baryon chemical potential is kept 
fixed at zero. In this case the Wroblewski factor is indeed seen to be a monotonic function of 



^ There the statistical model was fitted with an extra parameter 7s to account for possible chemical under- 
saturation of strangeness. At the SPS, 7s — 0.7 give the best agreement with Att data. 
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energy. The assumption of vanishing net baryon density is close to the prevaihng situation in 
e.g. p-p and e+-e~ colhsions. In Fig. 2 the results for extracted from the data in p-p, p-p 
and e"'"-e~ are also included |T^. The dashed line represents results with /i^ = and a radius 
of 1.2 fm. There are two important differences in the behavior of in elementary compared to 
heavy ion collisions. Firstly, the strangeness content is smaller by a factor of two. In elementary 
collisions particle multiplicities follow the values given by the canonical ensemble with radius 1.2 
fm whereas in A-A collisions there is a transition from canonical to grand canonical behavior. 
Secondly, there is no significant maximum in the behavior of A^ in elementary collisions due to 
the vanishingly small baryon density in the p-p and e"'"-e~ systems. 

The position of the maximum is further clarified in Fig. |^ which shows values of constant 
As in the T — plane. As expected A^ rises with increasing T for fixed /i^. Following the 
chemical freeze-out curve, shown as a thick full line in Fig. ^, one can see that A^ rises quickly 
from SIS to AGS energies, then reaches a maximum around /i^ ~ 500 MeV and T ^ 130 MeV. 
These freeze-out parameters correspond to 30 GeV lab energy. At higher incident energies 
the increase in T becomes negligible but fiB keeps on decreasing and as a consequence A^ also 
decreases. 

The importance of finite baryon density on the behavior of A^ is demonstrated in Fig. 4 
showing separately the contributions to (ss) coming from strange baryons, from strange mesons 
and from hidden strangeness, i.e., from hadrons like and 77. As can be seen in Fig. 4, the origin 
of the maximum in the Wroblewski ratio can be traced to the contribution of strange baryons. 
Even strange mesons exhibit a broad maximum. This is due to the presence of associated 
production of e.g. kaons together with hyperons. This channel dominates at low ^/s and loses 
importance at high incident energies. 

Next we study how the behavior of the Wroblewski ratio is reflected in specific particle 
yields. Of particular interest is the K yield which at high energies is responsible for almost 
80% of the total strangeness production while at lower energies this contribution decreases to 
50%, due to the associated production with hyperons. 

The energy dependence of the K~^/ 7r+ ratio measured at midrapidity is shown in Fig. 5. The 
model gives an excellent description of the data, showing a broad maximum at the same energy 
as the one seen in the Wroblewski factor. In general, of course, statistical-model calculations 
should be compared with 47r-integrated results since strangeness does not have to be conserved 
in a limited portion of phase space. A drop in this ratio for An yields has been reported from 
preliminary results of the NA49 collaboration at 158 AGeV @]. This decrease is, however, 
not reproduced by the statistical model without further modifications, e.g. by introducing an 
additional parameter 7^ ~ 0.7 [0. This point might be clearer when final data, also at other 
beam energies will become available. 

The appearance of the maximum in the relative strangeness contribution becomes also 
obvious when considering ratios which are more sensitive to the baryon chemical potential. 
Figure 6 shows the energy dependence of the A/vr^, S~/7r"*" and the fi/vr"'" ratios. As can be 
seen from the figure there is a very clear pronounced maximum especially in the A/vr"'' ratio. 
The relative enhancement of A is stronger then that of S~ or Q~. There is also a shift of the 
maximum to higher energies for particles with increasing strangeness quantum number. These 
differences appear as a consequence of the enhanced strangeness content of the particles which 
suppresses the dependence of the corresponding ratio on fiB- 



5 



4 Conclusions 



In conclusion, using the energy dependence of thermal parameters determined from an analysis 
of available data, we have shown that the statistical model description of relativistic heavy 
ion collisions predicts that the yields of strange to nonstrange particles reaches a well defined 
maximum near 30 GeV lab energy. It is demonstrated that this maximum is due to the specific 
shape of the freeze-out curve in the T — hb plane. In particular, a very steep decrease of 
the baryon chemical potential with increasing energy causes a corresponding decline of relative 
strangeness content in thermal systems created in heavy ion collisions above lab energies of 
30 GeV. The saturation in T, necessary for this result, might be connected to the fact that 
hadronic temperatures cannot exceed the critical temperature ^ 170 MeV for the phase 
transition to the QGP as found in solutions of QCD on the lattice. 

The maximum in the relative strangeness content is unique to heavy ion collisions: there is 
no equivalent behavior in the collisions of elementary particles. It has already been observed 
in the energy dependence of the j-n^ ratio, and is predicted to be even more pronounced 
in the kj-n^ ratio. The maxima in the 'ET j-n^ and r^^/vr^ ratios are expected to occur at a 
slightly higher beam energy and to be less pronounced. 

It is thus clear that in order to study collective effects related to strangeness production, 
one interesting area is the energy region around 30 A- GeV. Such a program is being planned 
for the near future at GSI. 
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Figure 1: Behavior of the freeze-out baryon chemical potential fiB (upper curve) and the 
temperature T (lower curve) as a function of energy. The temperature T as a function of beam 
energy is determined from the unified freezeout condition (E) / (N) = 1 GeV W^. 
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Figure 2: The Wroblewski ratio As (for definition see text) as a function of ^/s. The thick 
sohd hue has been calculated using the freeze-out values of the temperature and the baryon 
chemical potential. The dotted line has been calculated using /i^ = and only varying T. The 
dashed line has been calculated using a radius of 1.2 fm, keeping /xb=0 and taking the energy 
dependence of the temperature as determined previously. All calculations are performed using 
strangeness saturation 7,, = 1. 
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Figure 3: Lines of constant Wroblewski factor (for definition see text) in the T — plane 
(thin sohd hnes) together with the freeze-out curve (thick sohd hne) [|T^. 
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Figure 4: Contributions to the Wroblewski factor (for definition see text) from strange baryons 
(dotted line), strange mesons (dashed hne) and mesons with hidden strangeness (dash-dotted 
hne). The sum of all contributions is given by the full line 
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Figure 5: ji^^ ratio obtained around midrapidity as a function of from the various 

experiments. For the references for all data points see 0, . The full line shows the results of 
the statistical model in complete equilibrium. The value at RHIC was estimated using results 
from the STAR collaboration on the j-n^ and / ratios, assuming tt^/tt^ = 1.007. 
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Figure 6: Prediction for the A/tt^ (note the factor 5), S /vr^ and the ^7 /tt^ ratios as a 
function of ^/s. For compilation of data see 
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